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Abstract. We tabulate the abscissae and associated weights for numerical 
integration of integrals with either the singular weight function (— loga;) m 
for exponents m = 1,2 or 3, or the symmetric weight function cos(7ra/2). 
Standard brute force arithmetics generates explicit pairs of these values for up 
to 128 nodes. 



1. Methodology 



The paper provides abscissae Xi and weights Wi for Gaussian integration with a 
power of a logarithm in the integral kernel on one hand, 

,i n 

(1) / fixji-iogxrdx^y^wifixi). 

J i=l 

or with a cosine in the integral kernel on the other, 

(2) / f(x)cos(irx/2)dx&^2wif(xi). 

i=i 

The Wi and x, are computed with the standard theory from roots of a system of 
orthogonal polynomials p n with norm [6, 12, 14] 

(3) (/,<?) = f f(x)g(x)(-logx) m dx, 

Jo 

and 

/i 
7TX 
f(x)g(x) cos— da;, 

respectively. A set of orthogonal (monic) polynomials p n (x) is bootstrapped from 

(5) P-i(x) = 0; po(x) = 1; p n +i(x) = {x- a n )p n {x) - b n p n -i(x). 

[Dependence of polynomials and coefficients a and b on the parameter m in the 
case (1) is not written down explicitly here.] Multiplication of the recurrence with 
Pn or p n -i and using the requirement of orthogonality proposes to calculate the 
coefficients and polynomials recursively with 

(6) a n = -; 

\Pn,Pn) 
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It\ u n u \ x PniPn— l) , „■> 

(7) b = 0; b n = (n > 0). 

\Pn-l,Pn-l) 

Remark 1. In cases like (2) where the weight in the integral is an even function 
and the integral limits are symmetric, all a n are zero. 

The standard further steps are 

• normalization of the polynomials such that their norm is unity, 

fo\ * I \ — Pn\ x ) 

(8) P n (x) 



(9) 



V(Pn,Pn)' 

• computation of all zeros Xi of pn{%) at some degree N . 

• computation of the weights Wi by 



i xN ]P*N P*N+l( X i)P*N '( x iY 

where [x 7V+1 ]p^ r+1 and [x N ]p* N are the leading coefficients of the two poly- 
nomials after normalization, and where the prime at p' denotes the deriva- 
tive with respect to x. 

We obviously add no new aspect to the established theory. The benefit is to those 
readers who need explicit abscissae-weight pairs and have no access to a multi- 
precision numeric library. 



2. Logarithmic Kernel 

The first part of the results extends tables that have been published in the 
literature for exponent m = 1, namely by Anderson for N up to 10 [2], by Danloy 
for N = 10 and N = 20 [5], and by King for N = 20 and N = 30 [9]. 

Remark 2. The variable substitution x = e~ y changes the format to 
(10) / f(x){-\ogx) m dx= / f(e-y)y m e-ydy 



which is alternatively evaluated with Gauss-Laguerre quadratures [1, (25.4.38)] [3, 
13]. 

Integrals of the form (3) are calculated for the polynomials that appear in the 
recurrence (5) term-by-term with the aid of the moments fx [7, 2.722], 

(11) ^n,m= / X n (-\0gx) m dx 



(n + l)™+i 
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The first polynomials p„ !m (x) look as follows: 



(12) 




= X - 


-1/4; 


(13) 


P2.1 


= x 2 


- 5/7 * + 252 ; 


(14) 


P3,l 


= x A 


3105 2 5751 4679 
2588 X ' 16175 X 258800' 


(15) 


Pi, 2 


= x - 


-1/8; 


(16) 


P2,2 


= x 2 


19 217 
37 X + 7992 ' 


(17) 


P3,2 


= x 3 


1632663 2 5619807 1568083 
1695176 ^ ' 26487125 ^ 242168000' 


(18) 


Pl,3 


= x - 


-1/16; 


(19) 


P2,3 


= x 2 


13 493 
35 X+ 45360' 


(20) 


P3,3 


= x 3 


129197997 ^ , 4147011999 ^ 19126701359 
166534960 X ' 32526359375 X 8326748000000 



P2.1 in particular has been written down earlier [11]. Two generic values are 

(21) p hm = x-2- 1 ~ m ; 

2 m +l _|_ Qrn-\-^ ^m-\-l^ — l—m _j_ £m-\-lQ— 1— m 

(22) p 2 , m = X + y n + 1 - ^ m+1 X H 3>n+l _ 4 m+l " 

The results are summarized in the ASCII hies log_iV_m in the ancillary directory, 
where TV covers the range 3 to 128 and m covers powers from 1 to 3. Each line 
contains a pair (xi,Wi). For improved readability, a blank line is inserted after 
each block of 5 nodes. The numbers have been stabilized to the 30 digits shown 
by cranking up the internal representation of numbers in a Maple program to 270 
digits. 

Remark 3. Related approximative cubatures where polynomials are not only multi- 
plied by also added to the logarithm in the kernel have also been discussed [8, 4, 10]. 

3. Cosine kernel 
The tools to assemble (2) start from repeated partial integration of [7, 3.761] 

(23) I * m ™^=En k J^^X^ H-)^mM2i-\-m), 
for non- negative integer m. The even moments are therefore 

/l m 
i x 2m cos( OT /2)dz = 2(2m)! ER* (2m _ 2fc) , (2/7r) 2fe+1 

4 _ / — m+i.— m. 1 , 16\ 
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The odd moments are zero because the cosine is an even function. The monic 
orthogonal polynomials start 

(25) po 

(26) V 2 

(27) p 3 

(28) Pi 

and have parities p^ n (x) = (—l) n p n (x). 

The results are summarized in the ASCII files cosineJV in the ancillary direc- 
tory, where N covers the range 3 to 128. The numbers have been stabilized to the 
30 digits shown by an internal representation of numbers in a Maple program with 
650 digits. 

Only the values with positive Xi or xi = arc tabulated; the duplicates of the 
nodes at the negative abscissae (with the same weights) are not added explicitly. 
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1; pi=x; 

2 1 8 

x- -1 + —; 

3 tt 4 - 48tt 2 + 384 



X 



(TT 2 - 8)7T 2 



4 tt 4 - 78tt 2 + 672 2 TT 6 - 1147T 4 + 17287T 2 - 6912 
X ~ 7r 2 (7r 2 - 10) X + 7T 4 (7r 2 - 10) ' 



